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ABSTRACT 

The Vehicle Assembly Building (VAB) was constructed in the mid-1960s to house the Saturn V moon 
rocket while it was being assembled. Designed to withstand hurricanes and tropical storms, the VAB has 
a foundation consisting of 30,000 cubic yards of concrete strengthened by 4,225 steel rods driven 160 feet 
into limestone bedrock. The goal of the VAB Sway Investigation, which began collecting data in April 
2010 and ended in November 2012, was to quantify the displacement or sway of the VAB as a function of 
wind loading. Unlike other tall buildings whose height is much greater than their average cross-section 
dimension, the VAB is roughly a box with a height L = 526 ft. In the case of skyscrapers, building sway is 
characterized by a damped oscillation that may be on the order of 0.1 Hz to 0.5 Hz. The VAB, on the 
other hand, is an over-damped harmonic oscillator characterized by a simple displacement in the direction 
of the wind load. This conclusion is based on the fact that no oscillatory behavior was ever observed in 
the entire VAB sway data set. Even though the investigation and data collection continued for over 2 
years, the wind speeds observed during that time did not exceed 50 kts. Therefore, the final data analysis 
consists of an extrapolation formula, which expresses the magnitude of the displacement d as a function 
of wind speed v [kts], height above the ground h [ft], and number of standard deviations n above the 
mean: 


d(n,v,h) = a e hna (hv / L) 2 [in] 


( 1 ) 


where a = 0.0001591 [inkts' 2 ], b = -j5 lnlO = 5.1487 [in 1 ], and er= 0.092 [in], are determined from the 
data analysis. Table 1 is an evaluation of Equation (1) for various wind speeds and for n = 0, 1 , and 2. 
The meaning of confidence level is discussed in the Discussion and Summary section. 


Table 1. Equation (1) Evaluated for Various Confidence Levels C at h = 300 ft 


v [kts] 

(10-min peak) 

0 a 

(C = 50%) 

la 

(C = 84%) 

2 a 

(C = 98%) 

10 

0.005 

0.008 

0.013 

20 

0.021 

0.033 

0.054 

30 

0.047 

0.075 

0.12 

40 

0.083 

0.13 

0.21 

50 

0.13 

0.21 

0.34 

60 

0.19 

0.30 

0.48 

70 

0.25 

0.41 

0.66 

80 

0.33 

0.53 

0.82 

90 

0.42 

0.67 

1.1 

100 

0.52 

0.83 

1.3 

110 

0.63 

1.0 

1.6 

120 

0.74 

1.2 

1.9 
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VAB SWAY INVESTIGATION 


1 BACKGROUND 

The Vehicle Assembly Building (originally named the Vertical Assembly Building) was used to assemble 
American manned launch vehicles from 1968 through 2011. Completed in 1966, the VAB was originally 
built to allow for the vertical assembly of the Saturn V rocket for the Apollo program. It was then used for 
housing Space Shuttle external fuel tanks and flight hardware, and was where Space Shuttle orbiters were 
mated with their solid rocket boosters and external fuel tanks. The VAB is 526 feet (160.3 m) tall, 716 
feet (218.2 m) long, and 518 feet (157.9 m) wide. It covers 8 acres and encloses 129,428,000 cubic feet 
(3,665,000 m 3 ) of space. Entry to the bays inside the building is made through the four largest doors in the 
world. Each door is 456 feet (139.0 m) high and takes 45 minutes to completely open or close. Another 
large door on the north entry that leads to the transfer aisle was widened by 40 feet (12.2 m) to 
accommodate the Shuttle orbiter. A slot was added at the center of this north entry to allow for passage of 
the orbiter’s vertical stabilizer. 

It is well known that buildings sway in response to wind, as well as to earthquakes, and the taller the 
structure, the more they can sway. Because Kennedy Space Center (KSC) is occasionally assaulted by 
tropical or hurricane force winds, data on the VAB’s sway, collected from instrumentation positioned as 
indicated in Figure 1, is essential to help characterize the VAB facilities so that they can be used by 
commercial customers in the future. 



Figure 1. Locations for Data Collection 

Data from the nearest wind towers (http://trmm.ksc.nasa.gov/ ) is used to 
correlate wind velocity and direction to position data from the telescope. 
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2 TEST PROCEDURE 

The data acquisition strategy for this study was implemented with a Meade RCX400 telescope, positioned 
on the VAB floor in High Bay 4 as shown in Figure 2, along with a target fixed to the ceiling (see 
Figure 3). The position of the target is recorded once per second with a computer (also shown in Figure 2) 
attached to a camera mounted at the end of the telescope. The telescope data is then compared to data 
from the wind towers in the vicinity of the VAB (Figure 1). Instrumentation is located on three elevations 
on the wind towers; data from the highest, at 54 ft, is used in this analysis. The primary goal of this 
project is to correlate the movement of the building as a result of wind only. Table 2 is a list of the 22 
best-correlated data points extracted from telescope-target position data and the wind tower dataset. 



Figure 2. Telescope Used for Data Collection 


Meade RCX400 telescope pointed toward the ceiling in High Bay 4 with a PC 
running LabView image processing software. 


2 
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Table 2. Correlation of Wind Data with VAB Displacement 

(The correlation between wind data and VAB displacement is also shown as 
circles in Figure 6.) 


Date 

Wind speed 
(10-min peak) 

Displacement 
d[ in] 

Average Wind 
Direction [deg] 

10/27/2012 

39 

0.40 

0 

09/06/2012 

30 

0.12 

260 

09/09/2012 

27 

0.12 

273 

08/28/2012 

32 

0.08 

226 

08/27/2012 

35 

0.19 

167 

08/08/2012 

34 

0.21 

0 

07/22/2012 

41 

0.23 

167 

07/18/2012 

32 

0.15 

250 

06/25/20 12-a 

37 

0.18 

167 

06/25/20 12-b 

38 

0.23 

250 

06/06/2012 

42 

0.18 

250 

05/28/2012 

40 

0.18 

178 

04/06/2012 

29 

0.18 

0 

03/31/2012 

44 

0.54 

300 

03/04/2012 

36 

0.38 

310 

02/22/2012 

26 

0.27 

330 

09/14/2011 

23 

0.17 

47 

09/05/201 1 

31 

0.13 

250 

08/28/2011 

45 

0.13 

214 

06/04/2010 

31 

0.09 

226 

05/17/2010 

28 

0.18 

0 

05/06/2010 

35 

0.13 

250 
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Figure 3. Telescope Target 


The telescope’s target is mounted on the ceiling above the telescope in High Bay 4. 
As shown in the lower left, the target consists of three bright LEDs, in a rotated 
L-shaped pattern. The dimensions of the LED pattern are 1 inch x 3 inch. 


4 
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3 DATA ANALYSIS 

An example of correlated wind tower and telescope data is shown in the graph of Figure 4. At 
approximately 6 PM local time on August 8, 2012, a strong wind event occurred, lasting for well over an 
hour. The peak displacements (the differences between telescope and target positions) were only !4 inch, 
but were very well correlated to the wind tower data. The peak wind, as shown in Figure 4, was only 
35 kts, but the duration and broad impulsive character combined to generate a set of conditions that 
allowed the building position data to track the wind data. It is believed that a doubling of wind speed 
(70 kts in this case) will produce a quadrupling of the building sway, or 1 inch. This has not yet been 
confirmed, however; tropical force (70 kt) winds have not occurred since the hurricanes of 2004, or about 
5 years before the deployment of the sway-tracking telescope. 

In addition to wind-induced displacement, movement resulting from other sources, such as sun-induced 
heating and temperature cycling, as well as crane and heavy equipment operations, has been recorded. It 
is often difficult to distinguish the uncorrelated sources of telescope displacement noise, such as wind 
blowing through the large high bay doors. Since the telescope-target optical measurement reports a 
relative position, this measurement represents the maximum possible building sway, because it is 
statistically unlikely that the telescope and building would move in a way that would cancel out the real 
displacement. This sets a lower limit on displacement (but not an upper limit). In other words, the actual 
building displacement is equal to or less than the value measured by the displacement magnitude (< d x 2 + 
dy) m in the raw data plots, such as that shown in Figure 4. 
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Figure 4. Data from August 8, 2012 Event 

This wind event produced a well-correlated building displacement as measured by 
the telescope-target system. The wind direction was from the north, i.e., 0 degrees. 


In most of the good correlations between VAB displacement and wind events, the wind events are short 
and usually associated with thunderstorms or other phenomena related to passing weather fronts. The 
August 28, 201 1 data shown in Figure 5 is a prime example of a passing convective thunderstorm. When 
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this occurs, it is common for the wind speed to increase dramatically, with an abrupt change in direction. 
(Directional change is usually characterized by a rotation of wind direction.) After the front (sometimes 
referred to a gust front or outflow boundary) passes, the wind speed and direction return to their previous 
synoptic state. The entire length of this event is on the order of one hour or less, with the peak wind state 
lasting only 1 0 to 15 minutes. 

VAB Sway Wind Tower: AUG 2011 



-1000 




Figure 5. Example Thunderstorm Event and Associated Wind Direction Change over the 

VAB 


The thunderstorm event occurred on August 28, 2011, 17:58 UTC. Note that the 
time scale on the plot corresponds to local time. 
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4 DISCUSSION AND SUMMARY 


Figure 6 is a graph of Equation (1) as a function of wind speed for various values of n , and for h = 300 
(green lines) and h = L (red lines). When the data (red circles) is viewed on a log-log plot, such as 
Figure 6, the scatter appears to be uniformly distributed about the mean (0-a) line, shown as the bold red 
line, a is calculated by first transforming the wind speed v and displacement building d to new variables 
w = log(v) and u = log(d). The transformed data {w, w(w)} in the coordinate system is then rotated by 

6 - tan" 1 2 , since the slope of the log-log plot is 2, corresponding to the initial model, which assumes 
that displacement is proportional to v 2 . The displacement data in the rotated coordinate system can be 
modeled as a random variable x, which is approximately uniformly distributed about x = 0: 



( 2 ) 


where the standard deviation a is computed in the standard way in this rotated coordinate system. In order 
to use a in the original coordinate system of {v, d{y)}, an inverse rotation is performed followed by an 
anti-log transformation, arriving at: 

d(n,v,h) = a e bn<r v 2 [in] , (3) 

where a and b are fit parameters. By multiplying Equation (3) by the height function, (/z/Z,) 2 , which 
accounts for the quadratic bending mode of the building. Equation (1) is obtained. 



Figure 6. Building Displacement Correlated with Wind Speed 

Figure 6 shows Equation (1) evaluated at the (red) top of the VAB (h = L = 526 ft) and 
(green) at h = 300 ft. The solid bold lines are the average response (C = 50% confidence 
level), the solid line is one standard deviation above the mean (C = 84% confidence level), 
and dotted lines are two standard deviations above the mean (C = 98% confidence level). 
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The confidence levels C in Table 1 are obtained by integrating Equation (2) from -oo to no: 

MT 

C (n) = J P(x)dx 



The confidence level C can then be interpreted in this context as the percent fraction of the distribution 
based on the telescope-target and wind tower data, such that a building displacement at a given 
confidence level (based on n) is less than or equal to t/of Equation (1), corresponding to the fraction of 
P(x) from -oo to no. 

For this to have a more physical explanation, we need to discuss the error sources that contribute to the 
distribution and the measured standard deviation cr. The raw data can be decomposed into two main 
sources: (1) displacement that is correlated to the wind speed and direction at the tower, and 
(2) everything else. The second category corresponds to the sources mentioned in the previous section. 
The data used in this analysis (shown by circles in Figure 6) are only from the first category, i.e., where 
displacement is proportional to the wind speed and direction at the tower. For the purpose of this analysis, 
everything else was discarded. So the next question might be. What sources control the standard deviation 
cr, especially since all uncorrelated data was discarded? Many possibilities must be considered. 

The first source that contributes to the probability distribution width is the likelihood that not all of the 
uncorrelated noise was removed from the data used in Figure 6. This is simply because the correlated 
portion of the displacement is riding on a background signal that is uncorrelated. So even though attempts 
have been made to remove that background noise, some of it will remain, or the attempt at removing (or 
filtering) leads to additional error. 

However, we believe that the largest contribution of data scatter, quantified by cr (probability distribution 
width), can be investigated by answering the question. What is the actual wind speed acting on the VAB? 
Even though the wind towers are part of a quality-controlled weather data network for KSC and Cape 
Canaveral Air Force Station, they measure the wind at a specific point in space. Three wind towers were 
used for this analysis, as shown in Figure 1. The distance from the VAB to the towers is in a range from 2 
km to 3 km. The wind tower measurement used in this analysis was taken at a height of 54 ft (the highest 
available sensor on the three towers used). 

Another source of data scatter is the width of the time interval of the wind acting on the VAB. If the time 
interval is smaller than the response time of the building, the displacement will be smaller than the 
average trend. This could explain the circles under the main trend line (shown as the bold solid red line) 
in Figure 6. The circles above the line could be explained by wind forces at 300 ft that were greater than 
those measured by the wind towers at 54 ft. This is exemplified in the August 28, 201 1 event shown in 
Figure 5, where the entire event lasted less than one hour and the peak winds were prevalent for less than 
15 minutes. In this case, the building response time was longer than the wind event duration, and 
corresponds to a data point below the main trend line. 

So the primary uncertainty is in determining actual wind speed across the surface of the VAB. 
Conceptually, this could be improved by wind towers at 0.5 km (not too close, else the building-induced 
wind circulation would perturb the measurement), and at a height of 300 ft. A SODAR wind measurement 
system near the VAB could, in theory, produce an ideal wind measurement dataset. On the other hand, the 
primary data source used to monitor wind speeds relevant to VAB operations should be the ones used to 
make the associated predictions of building displacement. If that primary source is wind towers near the 
VAB (within 2 km to 3 km), then the resulting data scatter shown in Figure 6 must be dealt with in an 
appropriate statistical manner for useful displacement predictions, as has been attempted in this report. 
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APPENDIX A. RAW DATA CORRESPONDING TO TABLE 2 

The following 22 plots correspond to the entries of Table 2. Each plot contains wind speed and direction 
measured at the three towers shown in Figure 1. The horizontal and vertical target positions are also 
plotted, in black and grey. In order to extract a building sway direction, we need to approximate a local 
background position and apply it to the specific event. Flat lines in the telescope data usually indicate the 
crane blocking the telescope. The relative displacement of the building can then be converted into an 
angle for comparison with the wind direction recorded at the tower. Lastly, the LED error is plotted in 
red. 

The target position is the average of the three LED positions, (xo^o), (*i,>>i), and (^ 2 ,^ 2 ) over a 5-minute 
period: 


1 N 

c * = c xo + — X (*o (») + (») + x, (»)) (A- 1 ) 

1 N 

Cy=Cy0+ w^} y ° ( ' n ^ + y ' + y ' ^ ’ ( A ’ 2) 

where and Cyo are constant offsets for the graph. 

The LED error E (red line) is computed as the maximum of the 1 sec error over a 5-minute period, so that 
N=300: 


do\={( x o-X\) 2 +(y 0 -yi) 2 1 2 

(A-3) 

d n =((*, -X 2 ) 2 +(y { -y 2 ) 2) f 2 

(A-4) 

d 2 o = ((*2 - * 0 ) 2 + (y 2 - y 0 ) 2 1 2 

(A-5) 

E = ( h , - 1) 2 + (d l2 - 3) 2 + ( d 20 -\0 v2 ) 2 f 2 

(A-6) 


where the LED separations are assumed to be 1, 3, and 10 1 2 inch. A high error (under 0.5 inch) might 
indicate telescope movement, possibly from wind blowing past it from the high bay doors. An error 
higher than 0.5 inch usually indicates a missing LED (i.e., the camera does not detect it for a number of 
possible reasons). In some cases the error disappears from the plot for many days or weeks at a time. This 
occurs when the LabView imaging software fails to detect one or two LEDs. Fortunately, the 
measurement of the relative position of the target is not severely affected and useful building sway 
information can still be extracted. 
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APPENDIX B. 


DEVELOPMENT OF THE SWAY MATH MODEL 


Development of Equation (1) begins with the following physical principles as starting points: the sway 
displacement of the building d is a quadratic function of the wind speed v and is separable from the 
quadratic dependence of the vertical height h of the sway measurement: 


d(v,h) = f(v)g(h) , 


(B-l) 


/(v) = orv 2 


(B-2) 


g(h) = 0h 2 


(B-3) 


where a is parameter dependent on the building structure. The parameter p is equal to L' 1 where L is the 
building height. 

The data from Table 2 is plotted in Figure B-l, along with Equation (B-l) with g(L ) = 1. The parameter a 
is chosen to fit the data. The key to the strategy leading to the sway model represented by Equation (1) is 
finding the standard deviation of the data, corresponding to Table 2. 



20 30 50 70 100 

Wind Speed [kts] 

Figure B-l. The data from Table 2 with Equation (B-l) withg(Z) = 1. 

The first step is to rotate the coordinate system so that the mean of the data equals zero. The following 
variable substitution will simplify this procedure (where “log” is base 10): 

u — log d — u 0 , (B-4) 

w = logv , (B-5) 

where u 0 = logar . 


23 


N ASA/TM— 20 1 3-2 1 6320 


Figure B-2 shows Table 2 data plotted using the transformed variables (w w). Note that both plots are 
practically identical, where only the axis labeling and plot range has changed, as well as vertical offset, 
wo- 



W 

Figure B-2. Data from Table 2 plotted using variable transformations, Equations (B-4) and (B-5). The 

solid line has a slope of 2 and is just u(w) = 2w . 


The vector formed by ( w u ) is now rotated by an angle 9= tan' 1 2, which has the effect of forcing the 
slope of the fit in Figure B-l to zero: 




K u j 


{ cos 6 

sin<9" 

(w) 

^-sin# 

cos#. 

u, 


Figure B-3 is Table 2 data after application of Equation (B-6). 


(B-6) 
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•3C % * 










L0 0 1 2 3 4 



w' 

Figure B-3. Data from Figure B-2 plotted after undergoing the rotation transformation of Equation (B-6). 
The bottom plot is a zoomed-in version of the top plot. The line fit from the previous plot, u{w) = 2w, is 

now m'(V) = 0. 



Figure B-4. The u' k sequence from Figure B-3 with same ordering shown in Table 2. 


25 


N AS A/TM— 20 13-216320 


The next step is to compute the mean and standard deviation of the data in Figure B-3. Before doing so, it 
is helpful to replot Figure B-3 as a sequence of data u k , dropping the w k , since the random variable of 

interest is u k . Figure B-4 is the u k sequence from Figure B-3, in ascending order, the same ordering 
shown in Table 2. 

Now it is a simple matter to calculate the mean // and standard deviation a of the sequence of N data 
points shown in Figure B-4: 


The mean must equal zero, and it is necessary to iteratively adjust a , where u 0 = loga , so that ju —> 0. 
The result is shown in Table B-l. 


N 



(B-7) 



(B-8) 
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Table B-l. Adjusted data plotted in Figure B-4, with a= 0.0001591. 


k 

K 

1 

0.0976349 

2 

- 0.0342895 

3 

0.0066372 

4 

- 0.1381 1 

5 

- 0.0049171 

6 

0.0257814 

7 

- 0.0292712 

8 

- 0.0160197 

9 

- 0.037004 

10 

0.000245142 

11 

- 0.08624 

12 

- 0.0672877 

13 

0.0576297 

14 

0.109065 

15 

0.118765 

16 

0.178798 

17 

0.13657 

18 

- 0.0314805 

19 

- 0.176244 

20 

- 0.102901 

21 

0.0712608 

22 

- 0.0786225 

v 

-1.905 x 10" 10 

a 

0.0922 


Table B-l and Figure B-4 depict a transformed coordinate system where the sway measurement “noise” is 
independent of wind speed and only dependent on the errors resulting from random unknowns such as the 
wind correction factor between wind tower height and VAB height. Figure B-3 hints that as the horizontal 
variable w increases, the running average mean may not be zero. However, without higher wind speed 
measurements than what was observed during the period of this investigation, it is not constructive to 
infer a wind speed dependence in Figure B-3. But this could be worth reinvestigating for future analysis 
refinement. 

To get back to the original coordinate system where building displacement is a function of wind velocity, 
all of the previous transformations must be reversed. The first step on the way back is to express the 
function in Figure B-4 as a vector composed of a random variable x with zero mean, and a value above 
the mean given by no r. This vector is then acted upon by a rotation matrix that is the inverse of the 
original rotation: 
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^cos 9 -sin£^ 


( x' 

\U ) 


v sin<9 cos 9 j 




r xcosd-nusmB^ 
v /7<rcos# + ;csin#, 


The two components can be combined, eliminating x : 

u = ncr cos 9 + w” tan 9 . 


The last step is to add u 0 and take the anti-log (and dropping the double primes): 


d(n v) = io" < tW+mW+ “<> 

j Qwrsec0+logvtan0+loga 

= a lO^v 2 


= ae bna v 2 


Equation (B-ll) is now Equation (B-2). Note that b = y[5 In 10, since sec 9 = sec (tar 
Including the height dependence, Equation (B-3): 


d{n,v, h) = a e hn<T v 2 


r h b 2 




yields the desired result, Equation (1). 


(B-9) 


(B-10) 


(B-ll) 

'2)=V5. 

(B-12) 
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